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Abstract. The observation of a gamma-ray line in the cosmic-ray fluxes would be a smoking- 
gun signature for dark matter annihilation or decay in the Universe. We present an improved 
search for such signatures in the data of the Fermi Large Area Telescope (LAT), concentrating 
on energies between 20 and 300 GeV. Besides updating to 43 months of data, we use 
a new data-driven technique to select optimized target regions depending on the profile 
of the Galactic dark matter halo. In regions close to the Galactic center, we find a 4.6<r 
indication for a gamma-ray line at ~ 130 GeV. When taking into account the look- 
elsewhere effect the significance of the observed excess is 3.2<t. If interpreted in terms of 
dark matter particles annihilating into a photon pair, the observations imply a dark matter 
mass of m x = 129.8 ± 2.4 GeV and a partial annihilation cross-section of (av) xx -+~~ = 
(1.27 ± 0.32 toll) xl0~ 27 cm s when using the Einasto dark matter profile. The evidence 
for the signal is based on about 50 photons; it will take a few years of additional data to 
clarify its existence. 
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1 Introduction 

The existence of dark matter (DM) in the Universe is a by now well established fact, supported 
by numerous of observations that range from the time of matter / radiation equality until to- 
day [1]. A plethora of possible particle physics scenarios for DM exists, among which the 
currently leading hypothesis are weakly interacting massive particles (WIMPs) that are ther- 
mally produced in the early Universe [2-4]. WIMPs could be produced at particle colliders 
like the CERN LHC, they could be observed via WIMP-nucleon scattering in low-background 
experiments, and, due to their self-annihilation in the Universe, they could contribute to the 
cosmic-ray fluxes (see Ref. [5] for a recent review). 

Among the different particle species produced in WIMP annihilation, gamma rays are 
of particular importance, as they propagate essentially unperturbed and preserve spatial 
information about their sources. This information is exploited for efficient signal/background 
discrimination in searches for gamma-ray signatures from e.g. dwarf galaxies [6-14] or the 
Galactic center (GC) [15-18]. 

Besides spatial signatures, signatures in the energy spectrum of the gamma-ray flux play 
a major role in DM searches. It was noted long ago [19] that the two-body annihilation of 
DM into photons produces monochromatic gamma rays that could stand out of the otherwise 
continuous astrophysical fluxes as a clear smoking gun signature [20-26]. In most cases, 
the line signal is expected to be extremely faint, as its production is one-loop suppressed. 
However, it can be enhanced in many scenarios, like non-minimal variations of the MSSM [27, 
28], singlet DM [29], hidden U(l) DM [30, 31], effective DM scenarios [32], scenarios with 
strong annihilation into the higgs boson and a photon [33] or inert higgs doublet DM [34]; 
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these gamma-ray lines would be indeed in reach of current technology [35]. Additional 
prominent spectral features in the gamma-ray flux are generated by the emission of internal 
bremsstrahlung (IB) photons during DM annihilation [36-40]. 

Since more than three years, the Fermi Large Area Telescope (LAT), which is the main 
instrument on the Fermi Gamma-ray Space Telescope [41], measures the gamma-ray sky 
with an unprecedented precision, searching — amongst others — for DM signals [42]. Previous 
searches for spectral signatures from DM annihilation [43-47] were exclusively concentrating 
on gamma-ray lines (an earlier analysis using EGRET data can be found in Ref. [48]). A 
first dedicated search for IB signatures was recently presented in Ref. [49]; in this paper, a 
week indication for a spectral signature close to the GC was reported. 

The purpose of the present work is to update and refine the existing searches for gamma- 
ray lines from DM annihilation in Fermi LAT data. In particular, we are interested in a 
further investigation of the signature reported in Ref. [49]. Besides taking into account all 
data collected so far (43 months), and making use of the most up-to-date publicly available 
event selections, we employ a new data-driven algorithm to find target regions that maximize 
the signal-to-noise ratio (SNR) for different profiles of the Galactic DM halo. Since we are 
looking for very faint signatures, this is extremely important. An inefficient choice of the 
target region can easily hide a DM signal. We concentrate on line energies between 20 and 
300 GeV and use the data below 20 GeV for target region selection. 

The rest of this paper is organized as follows: In Section 2, we explain in detail the DM 
signal we are looking for, the treatment of the Fermi LAT data, the target region selection and 
the spectral analysis; in Section 3 we present our main results; a discussion of instrumental 
systematics, further tests and comparison with previous work can be found in Section 4; we 
conclude in Section 5. Event tables are provided in Appendix A. 



2 Data Analysis 

2.1 DM annihilation in the Galactic halo 

The gamma-ray flux from DM particles x annihilating into a photon pair inside the Galactic 
DM halo is (for a Majorana WIMP) given by 

^ K) = M^ M(£ _ £7) £ fp3m(r) , (2 . 1} 

where £ is the angle to the GC. Here, m x is the DM mass, {av) xx ^^ the partial annihila- 
tion cross-section for \X ~~ ^ 77) E~ = m x the gamma-ray line energy and Pdm( r ) the DM 
distribution as function of the Galactocentric distance r. The coordinate s > runs along 
the line-of-sight, and r(s,£) = y/ {tq — scos£) 2 + (ssin£) 2 , where tq = 8.5 kpc denotes the 
distance between Sun and GC. 

We consider five reference DM profiles, all normalized to /?dm( r o) = 0.4 GeV cm -3 at 
position of the Sun [50, 51]: (1) The Einasto profile, favored by the latest iV-body simula- 
tions [52-54], 

/ 2 r aE \ 

PdmM oc exp 5^ , (2.2) 

with cce = 0.17 and r s = 20 kpc. (2) The commonly used cored isothermal profile 

Pdm(r) oc 2 , (2.3) 

1 + (r/r s ) 
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with r s = 3.5 kpc (see e.g. Ref. [55] for observational arguments in favor of cored DM profiles). 
(3-5) The generalized Navarro-Frenk- White (NFW) profile 



(r/r s ) Q (1 + r/r s ) 

with r s = 20 kpc. Here, a parametrizes the profile's inner slope, and a = 1 reproduces 
the usual NFW profile [43, 56]. The presence of the super-massive black hole [57] or 
adiabatic contraction [58-61] can steepen the inner profile, which is here effectively taken 
into account by allowing that a > 1. To this end, we consider the three reference values 
a = (1.0, 1.15, 1.3), which are compatible with microlensing and dynamical observations [62]. 
Note that changing the slope of the inner profile in principle also affects the profile normal- 
ization. Since we are here mainly interested in the shape of the DM signal, however, we will 
keep Pdm(fo) = 0.4 GeV cm -3 fixed for simplicity. 

2.2 Data selection 

In the present analysis, we take into account 43 months of data (from 4 Aug 2008 to 8 Mar 
2012) with energies between 1 and 300 GeV. 1 We apply the zenith-angle cut 6 < 100° in order 
to avoid contamination with the earth albedo, as well as the recommended quality-filter cut 
DATA_QUAL==1. For comparison and cross-checks, we make use of both the SOURCE and 
ULTRACLEAN events selections (both Pass 7 Version 6). The former features an effective 
area that is relatively larger by about 12% at 100 GeV, the latter a lower contamination 
with residual charged cosmic rays (CRs). In the main analysis, front- and back-converted 
events are considered simultaneously, in order to minimize statistical errors. The selection of 
events as well as the calculation of exposure maps is performed using the 6 Oct 2011 version 
of ScienceTools v9r23pl. 2 For everything else we use our own software. 

2.3 Target regions & observed fluxes 

The optimal target region for gamma-ray line searches maximizes the SNR and depends 
on both, the morphology of the DM signal and the morphology of the background flux. 
For a particular Galactic DM profile, the former is given by Eq. (2.1); the latter has to 
be in principle determined from a modeling of gamma-ray diffuse and point sources at 20- 
300 GeV energies. The basic strategy that we follow here is to approximate the background 
morphology above 20 GeV by the spatial distribution of gamma rays measured between 1 
and 20 GeV. 

From events in the energy range 1-20 GeV, we produce a two-dimensional count map 
that covers an area of |fe| < 84° Galactic latitude and \£\ < 90° Galactic longitude, using 
a cylindrical equal-area projection with one square-degree pixel size. 3 For each pixel i, 
we derive the number of expected signal events \ii from Eq. (2.1), whereas the number of 
actually measured events is denoted by q. We need [ii only up to an overall normalization, 
since we leave (<7t> ) xx _> 77 and m x unspecified at this point. As long as the signal is only 
a weak perturbation, the SNR in each pixel can be estimated by TZi = Hi/yfci. Note that 
throughout the entire analysis, the angular resolution of the LAT — A9 ~ 0.2° for 20 to 
300 GeV energies [63] — is neglected. 



1 High level data of the LAT is available at http://fermi.gsfc.nasa.gov. 

2 See http://fermi.gsfc.nasa.gov/ssc/data/analysis for the standard chains. 

3 Because of this projection we exclude data with |6| > 84°. This does not affect the results. Note that we 
use the convention I = —180° . . . 180° throughout. 
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Figure 1. Left panel: The black lines show the target regions that are used in the present analysis in 
case of the SOURCE event class (the ULTRACLEAN regions are very similar) . From top to bottom, 
they are respectively optimized for the cored isothermal, the NFW (with a = 1), the Einasto and the 
contracted (with a = 1.15, 1.3) DM profiles. The colors indicate the signal-to- background ratio with 
arbitrary but common normalization; in Rcg2 to Rcg5 they arc respectively downscalcd by factors 
(1.6, 3.0, 4.3, 18.8) for better visibility. 

Right panel: From top to bottom, the panels show the 20-300 GeV gamma-ray (+ residual CR) 
spectra as observed in Regl to Reg5 with statistical error bars. The SOURCE and ULTRACLEAN 
events are shown in black and magenta, respectively. Dotted lines show power-laws with the indicated 
slopes; dashed lines show the EGBG + residual CRs. The vertical gray line indicates E = 129.0 GeV. 



-4- 



We define the optimal target region as the set of pixels T for which the integrated SNR 



is maximized. To find T , we use the approximate but efficient algorithm from Ref. [49]. 
The resulting target regions, optimized for our five reference DM profiles, are shown in the 
left panels of Fig. 1 by the black lines. The colors indicate the signal-to- background ratio 
/ii/cj, with an arbitrary but common normalization. In case of a cored isothermal profile 
(Regl), the target region is largest and reaches up to latitudes of |6| ~ 84°. The smallest 
region (Reg5) corresponds to a compressed profile with inner slope a = 1.3, and contains 
the central 2° x 2° degree of the GC only. In most cases, the regions are more extended 
south from the GC. This is a consequence of a slight north/south asymmetry in the observed 
diffuse gamma-ray flux (see e.g. Ref. [64]). Note that the regions are only optimal as long 
as signal contributions are small. In the presence of a potential signal these regions can be 
further optimized, which we leave for future work. 

We extract from the LAT data the gamma-ray flux measured in each of the five target 
regions. The corresponding energy spectra between 20 and 300 GeV are shown in the right 
panels of Fig. 1 for both, SOURCE (black) and ULTRACLEAN (magenta) event classes. The 
residual CR contamination of the SOURCE event selection is best visible in Regl as a sizeable 
difference between the SOURCE and ULTRACLEAN fluxes. This is further illustrated by the 
dashed lines, which show the expected flux of residual CRs plus the extragalactic gamma-ray 
background (EGBG) for comparison [65]. Remarkably, in Reg3 and Reg4 a pronounced bump 
at energies around 130 GeV (indicated by the vertical dotted line) can be easily recognized 
by eye; this spectral feature will turn out to be the best candidate for a gamma-ray line in 
the Fermi LAT data between 20 and 300 GeV. 

2.4 Spectral analysis 

In order to search for gamma-ray lines in Regl to Reg5, we perform a shape analysis of 
the energy spectra shown in Fig. 1 (though with much smaller energy bins). For a given 
gamma-ray line energy E^, this analysis is done in a small energy window that contains E^. 
The exact positions of the energy windows adopted during the main analysis are shown in 
Fig. 2. Since the energy windows follow the gamma-ray line energy, this method is known as 
"sliding energy window" technique [35, 43, 44, 48, 49]. 

We can conveniently parametrize the boundaries of the energy window as 



The sizes of the adopted energy windows vary between e ~ 1.6 at low energies, which is a 
few times wider than the LAT energy resolution, and e ~ 3.0 at high energies, in order to 
compensate for the lower number of events. The choice of the energy window size is somewhat 
arbitrary, but depends in principle on the uncertainties in the background curvature, the 
effective area and on the available statistics. We will discuss below how a change of the 
window size affects the results. 

Within the adopted energy windows, we fit the spectra from Fig. 1 with a simple three- 
parameter model, 



(2.5) 





(2.6) 




(2.7) 
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Energy windows used in the spectral analysis 




Figure 2. Energy windows that we use for our spectral line search. In red we indicate the window 
that enters the fit at E~ = 129.0 GcV. 



Background fluxes are here approximated by a single power law with a free spectral index 7 
and normalization ft, whereas the monochromatic DM signal has a free normalization S > 
while its position E~ remains fixed during the fit. Note that, after fixing the experimental 
conditions and the profile of the Galactic DM halo, the annihilation cross-section (cru) xx _ s . 77 
is related to S by a straightforward rescaling. 

The best-fit model parameters (5*bf, ftbf, 7bf) are obtained by maximizing the likelihood 
function C(S,ft,j) = HiP(si\ui), where P(s\v) = u s e~ u /si is the Poisson probability distri- 
bution function; here, Sj (i/j) denotes the number of measured (expected) events in energy 
bin i. In general, Vi is a function of the model parameters and calculated by multiplying 
the above three-parameter model with the exposure of the target region, and convolving the 
resulting function with the effective energy dispersion of the LAT. The details of this calcu- 
lation are discussed below in section 2.5. In our analysis we use energy bins that are much 
smaller than the energy resolution of the LAT, such that the analysis becomes independent 
of the actual binning and hence effectively "unbinned" . 4 

The significance of a line signal for a given value of £L is derived from the test statistic 

TS = -2\n^-, (2.8) 

•'-best 

where £best = £(<Sbf> ftbU 7bf) IS the likelihood of a fit with DM contribution, and £ nu n is 
the likelihood of a fit without DM signal (with S = fixed during the fit). In absence of a 
line-like signature in the data, the TS is expected to follow a 0.5Xfc =0 + 0.5x^ =1 distribution. 5 
Limits on, as well as statistical errors of, the annihilation cross-section (crt>) xx _j. 77 are de- 
rived using the profile likelihood method [66]. For instance, a one-sided 95% CL (two-sided 
68.2% CL) limit is obtained by increasing/decreasing S from S^f an d profiling over ft and 7 
until -21n(£/£ b cst) equals 2.71 (1.0). 

In the spectral analysis we scan over a large number of line energies E~ and over different 
target regions to find the maximal TS value. This reduces the statistical significance of any 
observed excess due to the well-known look-elsewhere effect. We attribute 6.8 trials over a 

4 We use 300 logarithmic bins per energy decade. 

5 This is a consequence of S > 0. The probability distribution function of xl=o i s S(TS). 
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X? =2 distribution to the quasi-continuous scan over line energies from EL = 20 to 300 GeV. 
This follows from a Monte Carlo analysis as well as from a subsampling analysis of the LAT 
data in regions without a DM signal; details can be found below in Section 4.3. Furthermore, 
we analysed ten different target regions, from which only the five most interesting ones 
are shown in this paper. 6 These target regions are partial subsets of each other, but we 
conservatively treat them as being statistically independent. However, we do not attribute 
trials to the scan over SOURCE and ULTRACLEAN event classes, as these are obviously 
strongly correlated. 

In summary, we find that the significance of a maximal TS value TS mSuX can in good 
approximation be derived from 10 x 6.8 = 68 trials over a xt=2 distribution. In practice, one 
has to solve 

CDF(xL 2 ;^max) 68 = CDF(xt =1 ;a 2 ) (2.9) 

for a. Here, CDF(x|;x) is the cumulative distribution function, which gives the probability 
to draw a value smaller than x from a xt distribution. 

2.5 Details on the treatment of instrument responses 

For a given infinitesimal sky region, the full expression for calculating the expected number 
of events Vi in a certain energy bin E~ . . . Ef reads 

rE t r r n ^ „ „ .,dT dj 

j=f,b 



V; 



dE J dE' j\e D{E\Ef,e,j)A{Ef,e,3)^^ . (2.10) 



Here, 9 is the impact angle of photon events with respect to the instrument axis, dT/dO 
the observational time under this impact angle, A(E',9,j) the effective area, E' the true 
energy, j = f,b denotes front- and back-converted events, E the reconstructed energy, and 
D{E\E' , 9, j) is the energy dispersion of the LAT. We define the effective exposure as 



x eff (E>) = j deJ^ME'Aj)^. (2.ii) 



j=f,b 

The expected distribution of the conditional parameters 9 and j for events with energy E' is 

P(B , ]F/) _ A(E',6,j)dT/d9 

P ^\ E )= Xeff (E>) ' (2 - 12) 
and we define the effective energy dispersion as 

D eff (E\E') = d9J2 D(E\E',9,j)P(6,j,E') . (2.13) 
h j=f,b 

With these definitions, we finally obtain that the number of expected events in energy bin 
i can be calculated as 



Vi 



dE J dE'D eff (E\E')X eff {E')^ . (2.14) 



The other target regions correspond to a = 1.05, 1.1, 1.2 and 1.4 as well as the Fermi Bubble template 
from Ref. [67] (see also Ref. [68]). 
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We use this expression, together with Eq. (2.7), to determine the number of expected events 
V{ for our likelihood analysis. 

The exposure X e ff(E') for the different target regions is derived using the ScienceTools, 
which automatically integrate over impact angles and front- and back-converted events ac- 
cording to the flight history of the instrument. We first generated lifetime cubes via gt It cube 
for the 43 months of data (with the options dcostheta=0 . 025, binz=l and zmax=100). Then 
we used gtexpcube2 to generate exposure maps for the instrument response functions (IRFs) 
P7SDURCE_V6 and P7ULTRACLEAN_V6. These maps are integrated over the target regions to ob- 
tain X e ff(E'); during this step, we weight the exposure maps with the profile of the expected 
dark matter signal to obtain correct signal normalizations. 

The parametrization of D(E\E' ,9,j) can be obtained from the ScienceTools. 7 To deter- 
mine P(9,j\E'), we directly use the distribution of events that enter our analysis. Integrating 
over the —90° < I < 90° region, we find that at energies of 100 to 300 GeV 55% of the events 
are front-, the rest back-converted. The distribution of impact angles 9 is derived from the 
-90° < t < 90° data above 1 GeV; it peaks at 9 ^ 40° and becomes negligible at > 70°. 
We approximate P(9,j\E') by adopting these distributions for all energies E' and for all ob- 
servational angles £1. This is a reasonable approximation for the following reasons: The ratio 
of front- and back-converted events is not strongly energy dependent, above 1 GeV it varies 
at most by 7%; since the energy resolution for front- and back-converted events differs above 
1 GeV at most by 20% [63], this weak energy dependence is negligible for our discussion. By 
inspecting the lifetime cube for different sky regions we find that, for the long time interval 
we consider, all regions of the sky are observed under a broad distribution of values with 
similar shape. However, since the maximal variation of the energy resolution with the impact 
angle is about 28% [63], one might consider the above approximation on the distribution 
as problematic. We estimate the impact of this approximation by adopting instead the 9 
distribution of events within a radius of 10° around the galactic center. We find that the 
FWHM and 68% containment width of the energy dispersion changes by less than 1%, which 
is negligible for our purposes. 

Note that the approximations for P(9,j\E') affect the energy dispersion only, i.e. the 
exact shape of the line signal, and not the calculation of the effective area or the power-law 
fits to the background, for which energy dispersion can be neglected anyway. For a region 
with 10° radius around the galactic center, we compared results obtained with our routines 
with the results that can be obtained by using gtrspgen (which is part of ScienceTools and 
takes fully into account effects of energy dispersion, but does not allow the use of arbitrary 
target regions) and XSPEC [69], and found excellent agreement. 

3 Main Results 

For each of the spectra shown in Fig. 1 we perform a search for gamma-ray lines in the 
range E~ = 20-300 GeV as described above. The resulting TS values as function of the 
gamma-ray line energy E^ are shown in the left and right panels of Fig. 3 for the SOURCE 
and ULTRACLEAN event classes, respectively. In regions Reg2, Reg3 and Reg4, we find TS 
values that are surprisingly large, and which indicate a high likelihood for a gamma-ray line 

7 The functional form of the energy dispersion can be found e.g. on 

http : 1/1 ermi . gsf c .nasa. gov/ ssc/data/analysis/documentation/Cicerone/Cicerone_LAT_IRFs ; the 

files in refdata/f ermi/caldb/CALDB/data/glast/lat/bcf /edisp/ in ScienceTools v9r23pl contain the 
required parameters. 
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Figure 3. TS value as function of the line energy E 7 , obtained by analysing the energy spectra from 
the different target regions in Fig. 1. Left and right panels show the results for the SOURCE and 
ULTRACLEAN event classes, respectively. The inset shows a zoom into the relevant region. The 
horizontal gray dotted lines show respectively from bottom to top the la to 3cr levels after correcting 
for trials (without trial correction the significance is given by \jTSa). In the right panel, the gray 
crosses show the TS values that we obtain when instead adopting the target region and energy windows 
from Refs. [45, 46] with 43 months of data. 



at Ej ~ 130 GeV. The largest TS value is obtained in case of the SOURCE events in Reg4 
and reads TS = 21.4 (corresponding to 4.6a before trial correction). Taking into account 
the look-elsewhere effect as discussed above, the trial corrected statistical significance for the 
presence of a line signal in the LAT data is 3.2a. 

The fits that yield the highest significance for a line contribution are shown in Fig. 4 
for the regions Reg2, Reg3 and Reg4, and for SOURCE and ULTRACLEAN events. In 
the upper sub-panels, we plot the LAT data with statistical error bars, as well as the total 
predicted counts from the best-fit models with (red bars) and without a gamma-ray line 
contribution (green bars). The blue dotted line shows the line flux component alone (before 
averaging over the energy bins). Note that, in order to improve the readability of the plots 
and to calculate the indicated p- values and the reduced Xred' we rebinned the data to five 
times fewer bins than actually used in the spectral fits. 8 The lower sub-panel shows the count 
residuals after subtracting the model with line. In most of the regions, the spectral signature 
that is responsible for the large TS values can be easily recognized by eye. The number of 
signal events ranges between 46 and 88, the statistical significance between 2.8<7 and 4.6a; 
the p-values and residual plots confirm that the fits to the data are reasonable and do not 
exhibit systematic discrepancies at low or high energies. 

If we interpret the observed signature as being due to DM annihilation into a photon 
pair via XX ~~ 77; we can constrain the DM mass m x (which then just equals the line energy, 
m x = -E 7 ) and the partial annihilation cross-section {av) xx ^^. The corresponding values for 
processes like XX ~ * jZ, 7/1 follow from a straightforward rescaling [44] . The inset of Fig. 3 
shows a zoom into the most interesting region of the TS plot. 9 From there, one can read off 
the DM mass that best fits the data together with its error bars. From the region with the 
largest TS value, Reg4 SOURCE class, we obtain m x = 129.8 ± 2.4+J 3 GeV. The indicated 

8 In Fig. 4, we omitted incomplete bins at the right end of the energy window. When calculating Xrodi we 
use the c-statistic Y^i 2(m* — c + 2c^ log(cj//ij). 

9 To generate the inset, we did not use sliding energy windows but kept the position of the energy window 
fixed at the position that corresponds to the with the largest TS. 
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Figure 4. Upper sub-panels: the measured events with statistical errors are plotted in black. The 
horizontal bars show the best-fit models with (red) and without DM (green), the Wwe dotted line 
indicates the corresponding line flux component alone. In the lower sub-panel we show residuals 
after subtracting the model with line contribution. Note that we rcbinned the data to fewer bins 
after performing the fits in order to produce the plots and calculate the p-value and the reduced 
xl = x 2 /dof. The counts are listed in Tabs. 1, 2 and 3. 
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Figure 5. Best-fit values for the annihilation cross-section into a photon pair, as obtained for different 
DM halo profiles in the different target regions of Fig. 1, together with their 68.2% CL and 95.5% CL 
errors. Previous upper 95% CL limits are shown for comparison [45, 46]. All values are derived 
assuming that m x = 129.0 GeV. 
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Figure 6. 68.2% (dotted), 95.5% (dashed) and 99.7% CL (solid) band of the annihilation cross-section 
(tTu) xx _i. 77 obtained when using a circular target region of radius 10° centered at the indicated values of 
l c and b c . The left (right) panel shows a latitudinal (longitudinal) scan. We assumed m x = 129.0 GeV 
during the fits. Note that the abscissa scales are different in both plots. 

errors are respectively statistical and systematical, the latter being due to uncertainties in 
the overall energy calibration of the LAT, AE/E =t\^% [70]. 

In Fig. 5 we show central values and 68.2% CL and 95.5% CL errors for the annihilation 
cross-section (crf) xx _j. 77 as derived from Regl to Reg5, assuming different DM profiles and 
m x = 129.0 GeV. We stress that the DM profiles only affect how signal normalizations 
from the spectral fits translate into annihilation cross-sections, but not the actual spectral 
analysis itself. In most of the regions, non-zero values of the annihilation cross-section are 
preferred at more than 95.5% CL, in agreement with Fig. 3. In case of the cored isothermal 
or the contracted DM profiles, we find some tension between the annihilation cross-sections 
obtained from different target regions. In case of the Einasto and NFW profiles, however, the 
values are mutually consistent. Using SOURCE class events and the spectrum from Reg4, 
we find best-fit annihilation cross-sections of (cu)^ x _j. 77 = (1.27±0.32±g;^|) x 1(T 27 cm" s 
in case of the Einasto profile, and of (<ti/) xx _j. 77 = (2.27 ± 0.57±g;|?) x 10~ 27 cm 3 s 1 in case 
of the NFW profile. The systematic uncertainties are here derived from the effective area 
(about 10% [70]) and from the energy calibration. 

In order to test the locality of the observed signature, we extract the gamma-ray energy 
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Figure 7. Upper 95% CL limits on the annihilation cross-section for \X ~^ 77; f° r different DM 
halo profiles, as obtained for the SOURCE event class. We used the correspondingly optimized 
target regions from Fig. 1. The green crosses show previous limits from Refs. [45, 46], assuming an 
Einasto profile; for comparison, the gray crosses show the limits that we obtain when using the energy 
windows and target region from Refs. [45, 46] with 43 months of data. The gray dotted line indicates 
m x = 130 GcV. 



spectra from a large number of circular target regions with a radius of 10°. These regions are 
either centered along the Galactic disk with b c = 0° , or they are centered at l c = 0° from the 
Galactic north to south pole; t c and b c denote the central coordinates of the target regions. 10 
From each of these target regions, we derived the 68.2%, 95.5% and 99.7% CL upper- and 
lower-limits on the annihilation cross-section (av) xx ->jj, assuming = 129.0 GeV and using 
the SOURCE event class. As shown in Fig. 6, we find that at 99.7% CL non-zero values of 
the annihilation cross-section are only preferred close to the GC; the observed signature 
disappears when moving to larger values of \l c \ or \b c \. 

We show 95% CL upper limits on the annihilation cross-section (<7u) xx _ 5 . Ty in Fig. 7. For 
each DM profile, the limits are derived from the correspondingly optimized target region. For 
all considered DM masses, the presented limits are well below the 'thermal cross-section' of 
3 x 10 -26 cm 3 s~ 1 . As expected, the limits are weakest in case of the cored isothermal profile, 
and strongest for the contracted profile with a = 1.3. At a mass around 130 GeV, they show 
a pronounced bump that corresponds to the large TS values in Fig. 3. For the Einasto profile, 
we plot previous results from Ref. [45, 46] — based on two years of data without a dedicated 
target region optimization — as green crosses; except for values around 130 GeV, we improve 
these limits by factors of two and more. 

The signal evidence is based on 46-88 photons, cf. Fig. 4. Drawing strong conclusions — 
in particular on the spatial distribution or the exact spectral shape of the signature — is thus 
hindered by the low number of events collected so far. Furthermore, our findings could 
be affected by instrumental systematics, which are difficult to evaluate with public data 
and information only. However, we checked that obvious systematic uncertainties do not 
significantly affect our findings; this will be discussed in the next section. 

10 Note that for b c 7^ 0° these regions are circular with respect to the projection used in Fig. 1. 
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4 Discussion 



4.1 Instrumental systematics 

Energy reconstruction and line shape. The present analysis is based on the P7SOURCE_V6 
and P7ULTRACLEAN_V6 event classes. The corresponding publicly available energy infor- 
mation (CTBBestEnergy) derives from a classification tree analysis of the parametric correc- 
tion (PC) and shower profile (SP) reconstruction algorithms. With respect to the older event 
class P6DATACLEAN.V3, PC and SP were debiased; the LK (likelihood) reconstruction al- 
gorithm, which was found to introduce jagged peaks that could be mistaken for gamma-ray 
lines, was removed (see e.g. Ref. [46]). The data set underlying the present analysis is hence 
expected to be much better suited for line searches than previously available data. 

In previous line searches conducted by the Fermi collaboration (based on the SP energy 
alone) , the spectral shape of reconstructed gamma-ray lines was carefully studied by means of 
Monte Carlo (MC) simulations that were calibrated with high-energy electron beam tests [46]; 
good fits to the MC results were obtained by a superposition of three Gaussians, leading to 
a peaked line with extended tails. In our analysis, the line shape is derived from the IRF 
that ships with the ScienceTools; it has a different analytical form, is somewhat broader, 
and related systematic uncertainties are for us difficult to control. In order to estimate to 
what extent different spectral line shapes may affect our results, we took the line shape 
from Fig. 5.3 of Ref. [46], as well as a simple single Gaussian with a FWHM identical to 
the one derived from the public IRF (13.6%), and used them in the fits. We find that the 
TS value corresponding to the observed signatures in Reg4 (SOURCE) decreases by 0.7 and 
3.0, respectively. The respective best- fit cross-sections are lowered by 8% and 30%. A more 
careful modeling of the line shape than what we could do in the present work is hence indeed 
relevant as it obviously affects the results, but it is extremely unlikely that it would invalidate 
our findings. 

Effective area and residual cosmic rays. Uncertainties in the effective area of the LAT 
have a large point-to-point correlation. On scales corresponding to the width of a gamma- 
ray line, they are about 2% (and grow to 15% at larger scales) [70]. This is much smaller 
than the fractional contribution of the observed feature to the overall flux, and can be hence 
neglected to good approximation. In any case, a sharp line-like feature in the effective area 
(or residual contamination from the nearly isotropic CRs) would also affect regions of the 
sky away from the GC. In light of Fig. 6 and the below subsampling analysis, this option 
appears very unlikely. 

Event selection. We checked that the signature appears in both, front- and back- 
converted ULTRACLEAN events separately, with a higher significance in back-converted 
events. The signature grew over time, with TS = 2.4 (8.8, 16.9) when taking only into account 
only the first 53 (107, 134) weeks of data from Reg3 SOURCE class. Furthermore, we checked 
that our results remain practically unchanged when using the filter cut DATA_QUAL==1 
&& IN_SAA!=T && LAT_CONFIG==l && ABS(ROCK_ANGLE)<52 instead, which re- 
duces the number of events by 2% with respect to the adopted DATA_QUAL==1 cut. 

4.2 Energy window size 

As already discussed above, one critical component in the analysis is the choice of the energy 
window size. In Fig. 8, the black lines show the TS values obtained when adopting different 
window sizes e, assuming a gamma-ray line energy of = 129.0 GeV. The window borders 
are given by Eq. (2.6) above, and we show results for Reg2 to Reg4 and SOURCE and 
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Figure 8. The black line shows how the TS value changes as function of the adopted energy window 
size, e. The gray shaded areas are the 68.2% and 95.5% CL bands for the TS values obtained from 
a MC simulation. The energy windows borders are calculated according to Eq. (2.6). We assumed 
E y = 129.0 GeV when generating the plot. 

ULTRACLEAN events. The window size that we actually used during the main analysis 
is e ~ 2.58 and indicated by the blue star. As one can see from Fig. 8, an increase of the 
window size in reasonable ranges leaves the TS values approximately unchanged. On the 
other hand, a further decrease of the window size would also decreases the TS value. This is 
expected, because the energy range over which the background normalization is fixed by the 
data becomes smaller and less constraining in that case. 

To study in more detail how well the behaviour of the TS values follows the statistical 
expectations, we perform a simple MC analysis: For each of the six cases in Fig. 8, we derive 
the best-fit model with line from an energy window with size e = 6.64. From this model, 
we generate 400 mock data sets and refit them using smaller windows. The 68.2% CL and 
95.5% CL bands of the TS values obtained in this way are shown in Fig. 8 by the dark 
and light shaded gray areas. We find that the actual behaviour of the observed TS value as 
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Figure 9. Black: Probability distribution function of maximal TS values that we found when per- 
forming line searches in 25000 Monte Carlo samples without line signal as described in the text. Gray: 
results obtained from 21000 subsamples of the hemisphere away from the GC. The blue line shows 
the fitting function PDF max (i; x); see text for details. 

function of e follows very well the MC results. This supports our choice of a window size of 
e ~ 2.58 at = 129.0 GeV. Smaller windows would significantly decrease the statistical 
power of the analysis. 

Lastly, we checked that fixing the slope of the background power-law to 7 = 2.6 (2.5, 
2.7) results in TS values of 20.4 (19.4, 21.5) in case of Reg4 (SOURCE); this again has little 
impact on our findings. A value of ~ 2.6 is physically expected since the gamma-ray flux in 
Reg4 is dominated by diffuse photons from cosmic-ray proton collisions with the interstellar 
medium. 

4.3 Monte Carlo and subsampling analysis 

In order to estimate the statistical significance of the observed excess, we perform a simple 
Monte Carlo analysis of our line search. We generate 25000 mock spectra that follow a power- 
law with slope —2.3 (like in Regl SOURCE class), such that each mock data set has the same 
number of counts as observed in Reg3 SOURCE. In each of these spectra, we perform the 
above spectral analysis, and record the corresponding largest TS value. The distribution of 
these maximal TS values is shown by the black dots in Fig. 9 with statistical error bars. We 
fit this distribution with the probability distribution function (PDF) of the maximum from 
t trials over a xt=2 distribution, which reads 

PDF max (i; x) = ^CDF(xL 2 ; *)* , (4-1) 

leaving t as free parameters in the fit. Note that our t corresponds to the 'effective number 
of search regions' discussed in Ref. [71] (related to the number of upcrossings), and k = 2 is 
in our case the theoretical expectation for the tail of the distribution. Since we are mainly 
interested in this distribution tail, we only take TS values above 5 into account in the fit. 
However, as shown in Fig. 9 by the blue line, we find that t = 6.8 trials well reproduce the 
empirical result over the whole range. We checked that increasing/decreasing the energy 
window or sample size within reasonable values does not change these results. 

In order to confirm the above Monte Carlo results with real data, we perform a sub- 
sampling analysis of the LAT data in regions where no signal is expected. To this end, we 
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exclude all data from the hemisphere towards the GC. The events from the region away from 
the GC, \b\ < 84° and \£\ > 90°, are distributed in energy bins i, the number of events in 
each energy bin is called Sj. From these we generate 21000 test spectra; the expectation 
value for each energy bin is given by /ij = O.Olsj. The small prefactor ensures that statistical 
fluctuations within the original data do not destroy the tail distribution of the subsampling 
analysis. Like above, we record for each of the test spectra the largest TS value found. The 
results are shown in Fig. 9 by the gray bands (which indicate the statistical error bars). They 
are in excellent agreement with the Monte Carlo results. We checked that this agreement 
does not strongly depend on the energy window size or on the above prefactor; it holds in 
particular also in the case where the number of events per test spectra equals the number of 
events in Reg3. 

We emphasize that line-like artefacts in the effective area of the LAT or contamination 
with isotropic cosmic rays would most likely have shown up as disagreement between our 
Monte Carlo and subsampling results, since they should also affect the fluxes observed towards 
the anti-galactic center. 

4.4 Comparison with previous results 

The strongest previously presented upper limits on gamma-ray lines at an energy of = 
130 GeV read (av) xx ^ < 1.4 (1.0, 2.4) x 10~ 27 cm 3 s" 1 for an NFW (Einasto, cored 
isothermal) DM profile (at 95% CL) [45, 46], and are based on two years of data and an 
analysis of a large sky region. In case of the Einasto profile, these limits are marginally 
consistent with the different best-fit cross-section that we show in Fig. 5. In case of the NFW 
(cored isothermal) profile, however, there is a two sigma tension between previous limits and 
the best-fit value from Reg4. This could be read as a weak preference for the Einasto (or 
for a compressed) DM profile. However, it could also indicate systematics related e.g. to the 
shape of the gamma-ray line that biases the obtained cross-sections to high, or it could be 
due to statistical fluctuations in the different underlying data sets. In any case, these limits 
would set interesting and complementary constraints on the allowed DM profiles if the 130 
GeV signature is due to DM annihilation. 

In Refs. [45, 46], which is based on the SP energy only, the significance for a signal 
at Ej = 130 GeV was found to be TS < 1. As a test whether this conflicts with our 
results, we implemented the target region and energy windows from Refs. [45, 46] (without 
doing a point-source subtraction), and repeated their scan over the discrete line energies 
£ 7 = 30, 40, . . . 200 GeV, using 43 months of ULTRACLEAN data. In this way, we find the 
TS values (upper limits) that are shown by the gray crosses in the right panel of Fig. 3 (Fig. 7). 
Indeed, no strong indication for a signature at 130 GeV appears, whereas the resulting limits 
are in good agreement with — and slightly stronger than — the two- years results from [45, 46]. 
This illustrates that a continuous scan over of gamma-ray line energies, a careful selection 
of energy windows and an improved selection of target regions is extremely important when 
searching for very faint signatures. 

5 Conclusions 

We presented a refined search for gamma-ray lines with the Fermi Large Area Telescope. 
The main improvement with respect to previous analyses is the use of a new optimization 
algorithm that automatically selects target regions with the largest signal-to-noise ratio, 
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depending on the adopted Galactic dark matter profile and the LAT data. Besides, we 
updated previous results to 43 months of data and used the most recent public event selection. 

We found a 4.6cr indication for the presence of a gamma-ray line at ~ 130 GeV. 11 
Taking into account the look-elsewhere effect, the significance is 3.2a. By scanning through 
different regions of the sky, we find that the signature appears close to the Galactic center 
only. If interpreted in terms of dark matter annihilation into photon pairs, XX ~~ ^ 77? the 
observations constrain the dark matter mass to m x = 129.8 ± 2.4^ 7 3 GeV. The annihilation 
cross-sections derived from different target regions are consistent with each other in the 
case of the Einasto or NFW dark matter profiles. For these, we obtain best-fit annihilation 
cross-sections of (av) xx ^^ = (1.27 ± 0.32±o;2s) x 10~ 27 cm 3 s" 1 and {av) xx ^ yy = (2.27 ± 
0.57^Q5^) x 10 -27 cm 3 s — 1 , respectively. Assuming a thermal relic, this corresponds to a 
branching ratio into 77 final states of 4-8%, which is much larger than what the generic one- 
loop suppression suggests. We examined our statistical method with Monte Carlo methods 
and with a subsampling analysis of the LAT data in the \£\ > 90° hemisphere, and find 
good agreement with theoretical expectations. We performed several checks to exclude that 
obvious systematics of the LAT could invalidate our findings. In particular the fact that the 
signal peaks around the Galactic center makes a purely instrumental cause very unlikely. 

The observation of a gamma-ray line from dark matter annihilation in the data of the 
Fermi LAT is a fascinating possibility. However, a few caveats are in order. First, the 
presented analysis is based on publicly available data and information only, hence we cannot 
take into account all possible instrumental effects. Second, the signal evidence is right now 
based on around 50 photons; it will require a few years of more data to settle its existence on 
statistical grounds. Third, due to the low number of events it is difficult to study in detail 
the spatial and spectral characteristics of the signature, which might still leave room for some 
exotic astrophysical explanation. 

In any case, we showed that a careful selection of target regions is extremely important 
when searching for gamma-ray lines. We hope that our work motivates a refined search for 
such features in the Fermi LAT data, including all instrumental effects, an improved event 
selection and a target region optimization as proposed in this paper. 
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European 1231 Union FP7 ITN INVISIBLES (Marie Curie Actions,PITN-GA-2011-289442). 
This work makes use of SciPy [72], PyFITS 12 , PyMinuit 13 and IPython [73]. 

Note added. After submission of this manuscript, the two-years line search of the LAT 
collaboration was released [74]. The results were already presented on different conferences 
during last year [45], details can be found in Ref. [46]. We comment on them where appro- 
priate. 

11 We reported this signature before in Ref. [49] in context of internal Bremsstrahlung signals from DM. 

http: //www . stsci . edu/resources/ sof tware_hardware/pyf its 
13 http: //code .google . com/p/pymiimit 
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E bins 


# Counts 




Flux 


[VJC V J 






[ph/cm /s/ sr] 


79.3-82.4 


64 


0.34 


1.37e-10 


82.4-85.6 


84 


0.42 


1.73c-10 


85.6-89.0 


71 


0.53 


1.47e-10 


89.0-92.5 


66 


0.68 


1.27c-10 


92.5-96.1 


63 


0.88 


1.17c-10 


96.1-99.8 


68 


1.17 


1.22e-10 


99.8-103.7 


68 


1.61 


1.18e-10 


103.7-107.8 


65 


2.31 


1.09e-10 


107.8-112.0 


55 


3.54 


8.88e-ll 


112.0-116.4 


72 


5.81 


1.12e-10 


116.4-120.9 


50 


9.68 


7.53e-ll 


120.9-125.7 


46 


15.18 


6.7c-ll 


125.7-130.6 


67 


18.89 


9.43e-ll 


130.6-135.7 


51 


13.76 


6.94e-ll 


i 9tr 7 1/11 A 

loo. /— 141.0 


36 


7.07 


4.74e-ll 


141.0-146.5 


35 


3.04 


4.45c-ll 


146.5-152.2 


34 


1.30 


4.18c-ll 


152.2-158.2 


34 


0.61 


4.04e-ll 


158.2-164.4 


28 


0.34 


3.22c-ll 


164.4-170.8 


23 


0.20 


2.56e-ll 


170.8-177.5 


26 


0.12 


2.79e-ll 


177.5-184.4 


22 


0.07 


2.29e-ll 


184.4-191.6 


19 


0.04 


1.91e-H 


191.6-199.1 


11 


0.03 


1.36e-ll 



E bins #Counts #Signal Flux 
[GeV] [ph/cm 2 /s/sr] 



81.1-84.3 


47 


0.16 


9.39e-ll 


84.3-87.6 


51 


0.20 


9.83e-ll 


87.6-91.0 


IS 


0.25 


8.94e-ll 


91.0-94.6 


11 


0.33 


7.91c-ll 


94.6-98.3 


41 


0.43 


7.12c-ll 


98.3-102.1 


39 


0.58 


6.55e-ll 


102.1-106.1 


12 


0.81 


6.8e-ll 


106.1-110.3 


31 


1.16 


4.84e-ll 


110.3-114.6 


19 


1.78 


7.37e-ll 


114.6-119.1 


36 


2.94 


5.23e-ll 


119.1-123.7 


20 


4.96 


2.81e-ll 


123.7-128.6 


36 


7.99 


4.9e-ll 


128.6-133.6 


11 


10.57 


5.39e-ll 


133.6-138.8 


25 


8.16 


3.18c-ll 


138.8-144.3 


24 


4.25 


2.95e-ll 


144.3-149.9 


22 


1.82 


2.62e-ll 


149.9-155.8 


21 


0.76 


2.41e-ll 


155.8-161.9 


25 


0.35 


2.77c-ll 


161.9-168.2 


13 


0.18 


1.39e-ll 


168.2-174.8 


15 


0.10 


1.55e-ll 


174.8-181.6 


16 


0.06 


1.59e-ll 


181.6-188.7 


14 


0.04 


1.35e-ll 


188.7-196.1 


11 


0.02 


1.03e-ll 


196.1-203.8 


10 


0.01 


9.02e-12 


203.8-211.7 


6 


0.01 


5.24e-12 



Table 1. From left to right, the tables show the energy bins, the number of observed counts, the 
best-fit signal counts, and the observed fluxes that appear in Fig. 4. The left table corresponds to 
Reg2 SOURCE class, the right table to Reg2 ULTRACLEAN class. Although in our main analysis 
we use five times smaller energy bins, the listed numbers are enough to reproduce our results with 
good accuracy. 



E bins 


# Counts 


#Signal 


Flux 


E bins 


#Counts 


#Signal 


Flux 


[GeV] 






[ph/ cm 2 /s/sr] 


[GeV] 






[ph/cm 2 /s/sr] 


80.5-83.7 


31 


0.26 


2.78e-10 


80.5-83.7 


21 


0.18 


2.03e-10 


83.7-86.9 


33 


0.33 


2.86e-10 


83.7-86.9 


23 


0.23 


2.14e-10 


86.9-90.3 


25 


0.41 


2.09e-10 


86.9-90.3 


21 


0.29 


1.89e-10 


90.3-93.9 


21 


0.53 


1.69e-10 


90.3-93.9 


IS 


0.37 


1.56e-10 


93.9-97.6 


32 


0.69 


2.49e-10 


93.9-97.6 


21 


0.49 


1.76c-10 


97.6-101.4 


24 


0.91 


1.8c-10 


97.6-101.4 


16 


0.66 


1.3e-10 


101.4-105.3 


31 


1.26 


2.25e-10 


101.4-105.3 


24 


0.92 


1.88e-10 


105.3-109.5 


17 


1.80 


1.19e-10 


105.3-109.5 


14 


1.34 


1.06e-10 


109.5-113.7 


26 


2.76 


1.76c-10 


109.5-113.7 


21 


2.06 


1.53e-10 


113.7-118.2 


32 


4.53 


2.09e-10 


113.7-118.2 


20 


3.41 


1.4e-10 


118.2-122.8 


21 


7.56 


1.33e-10 


118.2-122.8 


14 


5.76 


9.49e-ll 


122.8-127.6 


23 


11.87 


1.41e-10 


122.8-127.6 


17 


9.21 


l.llc-10 


127.6-132.6 


32 


14.82 


1.89e-10 


127.6-132.6 


26 


11.80 


1.65e-10 


132.6-137.8 


21 


10.83 


1.2e-10 


132.6-137.8 


15 


8.72 


9.2c-ll 


137.8-143.2 


22 


5.58 


1.21c-10 


137.8-143.2 


17 


4.45 


l.Olc-10 


143.2-148.8 


7 


2.40 


3.74e-ll 


143.2-148.8 


6 


1.89 


3.44e-ll 


148.8-154.6 


12 


1.03 


6.19c-ll 


148.8-154.6 


9 


0.79 


4.98e-ll 


154.6-160.6 


9 


0.48 


4.49e-ll 


154.6-160.6 


5 


0.36 


2.67e-ll 


160.6-166.9 


9 


0.27 


4.34e-ll 


160.6-166.9 


7 


0.19 


3.61e-ll 


166.9-173.4 


10 


0.16 


4.66e-ll 


166.9-173.4 


7 


0.11 


3.48e-ll 


173.4-180.2 


10 


0.09 


4.51c-ll 


173.4-180.2 


8 


0.07 


3.84e-ll 


180.2-187.3 


8 


0.06 


3.48e-ll 


180.2-187.3 


7 


0.04 


3.25e-ll 


187.3-194.6 


8 


0.03 


3.37e-ll 


187.3-194.6 


5 


0.02 


2.25e-ll 


194.6-202.2 


6 


0.02 


2.44e-ll 


194.6-202.2 


4 


0.01 


1.74e-ll 



Table 2. Same as Tab. 1, but for Reg3. 
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E bins 


# Counts 


II o * „_ „ 1 

#Signal 


Flux 


[GeV] 






[ph/ cm 2 /s/sr] 


80.5-83.7 


19 


0.22 


2.99e-10 


83.7-86.9 


26 


0.27 


3.95e-10 


86.9-90.3 


16 


0.34 


2.35e-10 


90.3-93.9 


13 


0.43 


1.84e-10 


93.9-97.6 


21 


0.56 


2.87c-10 


97.6-101.4 


13 


0.74 


1.72e-10 


101.4-105.3 


18 


1.01 


2.29c-10 


105.3-109.5 


11 


1.45 


1.35e-10 


109.5-113.7 


17 


2.21 


2.02e-10 


113.7-118.2 


23 


3.60 


2.64e-10 


118.2-122.8 


13 


6.02 


1.44e-10 


122.8-127.6 


15 


9.52 


1.61c-10 


127.6-132.6 


24 


12.29 


2.49e-10 


132.6-137.8 


13 


9.37 


1.3c-10 


137.8-143.2 


19 


4.93 


1.84e-10 


143.2-148.8 


5 


2.14 


4.68e-ll 


148.8-154.6 


7 


0.91 


6.34e-ll 


154.6-160.6 


7 


0.43 


6.12e-ll 


160.6-166.9 


5 


0.23 


4.23e-ll 


166.9-173.4 


6 


0.13 


4.91e-ll 


173.4-180.2 


7 


0.08 


5.53e-ll 


180.2-187.3 


5 


0.05 


3.82e-ll 


187.3-194.6 


3 


0.03 


2.22e-ll 


194.6-202.2 


4 


0.02 


2.86e-ll 


202.2-210.1 


2 


0.01 


1.38e-ll 



Table 3. Same as 



E bins 
[GeV] 


^Counts 


it o: „„ „ l 

#bignal 


Flux 

[ph/cm 2 /s/sr] 


80.5-83.7 


14 


0.15 


2.31e-10 


83.7-86.9 


20 


0.19 


3.18e-10 


86.9-90.3 


13 


0.25 


1.99e-10 


90.3-93.9 


10 


0.32 


1.48e-10 


93.9-97.6 


14 


0.42 


2e-10 


97.6-101.4 


7 


0.57 


9.68e-ll 


101.4-105.3 


15 


0.79 


2e-10 


105.3-109.5 


9 


1.14 


1.16o-10 


109.5-113.7 


13 


1.75 


1.61e-10 


113.7-118.2 


17 


2.88 


2.04e-10 


118.2-122.8 


9 


4.87 


1.04e-10 


122.8-127.6 


13 


7.81 


1.46e-10 


127.6-132.6 


19 


10.18 


2.06e-10 


132.6-137.8 


9 


7.68 


9.42e-ll 


137.8-143.2 


15 


3.96 


1.52e-10 


143.2-148.8 


4 


1.69 


3.92e-ll 


148.8-154.6 


6 


0.70 


5.67e-ll 


154.6-160.6 


4 


0.32 


3.65e-ll 


160.6-166.9 


4 


0.17 


3.52e-ll 


166.9-173.4 


5 


0.10 


4.24e-ll 


173.4-180.2 


5 


0.06 


4.09e-ll 


180.2-187.3 


4 


0.03 


3.170-11 


187.3-194.6 


3 


0.02 


2.3e-ll 


194.6-202.2 


3 


0.01 


2.22e-ll 


202.2-210.1 





0.01 






Tab. 1, but for Reg4. 
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